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Abstract 

The vortex free energy is studied in the random vortex world-surface model of the in- 
frared sector of SU (3) Yang-Mills theory. The free energy of a center vortex extending into 
two spatial directions, which is introduced into Yang-Mills configurations when acting with 
the 't Hooft loop operator, is verified to furnish an order parameter for the deconfinement 
phase transition. It is shown to exhibit a weak discontinuity at the critical temperature, 
corresponding to the weak first order character of the transition. 
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1 Introduction 



Yang-Mills theory in four space-time dimensions undergoes a deconfinement phase transition at 
non-zero temperatures. One order parameter used in practice to characterize this transition is 
the Polyakov loop, which can be interpreted as representing a single static quark, immersed in a 
thermal Yang- Mills background. Since the confining phase is characterized by the free energy Fg 
of such a quark being infinite, the Polyakov loop (P) = exp{—Fq) is expected to vanish below the 
critical temperature Tc , while it should be non-zero in the deconfined (high temperature) phase. 
The Polyakov loop is closely related to the so-called center symmetry, which is broken at high 
temperatures T > T^. Physically, this critical behavior is counterintuitive; usually, e.g. in most 
spin models, the high temperature phase is maximally symmetric and the symmetries are broken 
below the critical temperature. Furthermore, when formulated on a compact space with periodic 
boundary conditions (such as used in lattice calculations), Yang- Mills theory strictly speaking 
does not admit a single quark as a physical excitation, regardless of whether confinement is 
realized or not. 

Not least due to these issues, it is instructive to consider an alternative dual (dis-)order parameter 
to characterize the deconfinement phase transition, namely the free energy associated with a 
center vortex world-surface. Center vortices are defined by the property that they contribute a 
center phase to a Wilson loop when piercing an area spanned by the latter. The center vortex 
free energy is expected to behave in a manner which is dual to the behavior of Wilson loops: 

• The free energy of vortex world-surfaces extending into one space direction and the (Eu- 
clidean) time direction is expected to exhibit no area-law dependence on the vortex world- 
surface area at all temperatures T, if the corresponding dual (spatial) Wilson loops exhibit 
an area law for all T (this is the case for pure SU{N) Yang-Mills theory). 

• The free energy of vortex world-surfaces extending into two spatial directions is expected 
to behave in a manner dual to the behavior of corresponding (temporal) Wilson loops: No 
area-law dependence of the free energy on the vortex world-surface area at temperatures 
below the deconfinement transition temperature Tc, but behavior according to an area law 
above the critical temperature. 

The latter observation permits the definition of a (temperature-dependent) dual string tension 
a through the leading behavior of the excess free energy Fp in the presence of asymptotically 
large vortex world-surfaces P, 

Fp = -ln{Zp/Zo) a(r)Area(P), T > T, (1) 

which can serve as an order parameter for the deconfinement phase transition. Here, Zq denotes 
the conventional Yang-Mills partition function, whereas Zp denotes the partition function in 
the presence of the additional vortex P. 

Using the vortex free energy as a confinement order parameter was first suggested by 't Hooft 
who, initially working in the Schrodinger picture (using the Weyl gauge), defined a magnetic 
loop operator M ('if) in the continuum P via equal-time commutation relations with all spatial 
Wilson loops W(^'), 

M(^)w('r')MCr)^ = z('r,'r')w('r'). (2) 

Here, z{^,^') is an element of the center of the gauge group related to the integral linking 
number v{'^,'^') of the two loops ^ and in three space dimensions, 

z{^, = exp f ^')\ ( SU{N) gauge group ). (3) 
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The loop operator M of course is only defined implicitly by the commutation relation eq. (j2I). 
An explicit realization of this formal construction was recently derived and discussed in detail in 
ref. [21, cf. also Also in the dual language, the change of behavior of the vortex free energy 
across the deconfinement phase transition can be associated with a (magnetic) center symmetry 
^13], which is broken in the low-temperature phase. 

Subsequent to 't Hooft's suggestion, realizations of the 't Hooft loop operator within the frame- 
work of lattice gauge theory formulated in 3+1-dimensional (Euclidean) space-time were dis- 
cussed and schemes of measuring the free energies of the vortices thus introduced into the gauge 
configurations were devised El IZl E] • On this basis, sufficient conditions for confinement via 
vortex condensation were established IHl • While the construction given in admits the intro- 
duction of vortices of a large variety of shapes into gauge configurations (by imposing suitable 
boundary conditions on "vortex containers" within the lattice), in practice a particularly clean 
way to investigate the center vortex free energy is to consider specifically a vortex world-surface 
which occupies a complete two-dimensional plane in space-time^. On a compact space-time, 
such as used in lattice Yang-Mills theory, such a world-surface is closed by virtue of the periodic 
boundary conditions; furthermore, its introduction corresponds to adopting twisted boundary 
conditions in the directions orthogonal to the vortex JU] . This is reflected in the fact that config- 
urations with such a vortex are not connected to the conventional vacuum by the dynamics. The 
dynamics of Yang-Mills theory and also of the random vortex world-surface model considered in 
this work, cf. |1H I12j , can create at most pairs of vortices occupying complete two-dimensional 
(parallel) planes in space-time; a single such vortex is topologically stable. Also, its free energy 
is not contaminated by ancillary effects such as the ones which occur, e.g., when one consid- 
ers open vortex world-surfaces; these are bounded by center monopole world-lines, and their 
free energies thus contain additional contributions stemming from monopole self-energies and 
interactions, which need to be disentangled from the vortex free energy itself^. 
A number of measurements of vortex free energies have been carried out within lattice Yang- 
Mills theory, in a variety of guises. The excess free energy in the presence of twisted boundary 
conditions compared to periodic boundary conditions was evaluated in ^1 El ^1 El ITHl 
1191 1201 121j . Also the free energy of more varied vortex configurations, including open vortex 
world-surfaces (HI 1171 1221 123j and intersecting planar vortices was investigated. Partition 
functions in the presence of twists can furthermore be related to electric flux free energies via 
Z{N) Fourier transforms ^01; a thorough discussion of this together with detailed measurements 
can be found in jl8j . cf. also [241 125j for earlier work with this focus. In addition, in the case of 

^Note that the covariant setting given by the lattice gauge theory framework permits the definition of vortex 
world-surfaces extending in arbitrary directions in space-time; however, ones extending purely in two spatial 
directions are of course the relevant ones as far as defining an order parameter for the deconfinement phase 
transition is concerned. 

note is in order concerning the "center monopole" concept. Center monopoles in the strict sense, i.e., sources 
and sinks of magnetic flux of magnitude 2n/N, are unphysical objects, since they violate the Bianchi identity 
(continuity of magnetic flux modulo 2ti). In lattice Yang-Mills theory, when separating the gauge fields into 
center and coset components, one can identify locations at which magnetic flux carried by the center component 
spreads out and continues as part of the coset component |7|. These locations are often called "center monopoles", 
although strictly speaking they are not sources or sinks of magnetic flux and do not violate the Bianchi identity. 
On the other hand, in the vortex model discussed here, there is no analogue of the coset; center vortices (albeit 
implicitly endowed with a flnite thickness) are the only degrees of freedom, and open vortex world- surfaces, 
bounded by center monopoles, are unphysical. Nevertheless, further below, open vortex surface conflgurations 
will be used as intermediate objects in setting up the numerical calculation; these should not be viewed as any 
more than auxiliary mathematical constructs. The dual string tension measured will ultimately be derived from 
the free energy of a closed vortex world-surface. 
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a system with a first order deconfinement transition, as found, e.g., for SU{3) color, the vortex 
free energy in the deconfined phase is related to the Z{3) order-order interface tension P]; for an 
evaluation of this interface tension using a different method, cf. [211 ■ With auxiliary information 
about wetting properties |261 127j . the order-order interface tension can in turn be connected 
to the order-disorder interface tension at the critical temperature. Direct measurements of the 
order-disorder interface tension at the SU{3) deconfinement transition were reported in j271 1281 
E2I IHOl 1^ IH2I inni IM] ; for a recent comparison with data obtained using a new calculational 
method based on the vortex free energy, cf. |2J. The results obtained in the present work will 
be put in relation to the survey presented in |2j in section ID 

2 Vortex free energy in the random vortex world-surface model 

The main purpose of the work presented here is to compute the dual string tension d'{T) by 
evaluating the free energy of appropriate vortex world-surfaces within the SU{3) random vortex 
world-surface model. This model was first defined and studied for the gauge group SU (2) in 
refs. ^2 ESI ^'^d later extended to the SU{3) case in refs. It describes the infrared 

sector of Yang-Mills theory as an ensemble of random (thick) center vortex world-surfaces. On 
a space-time lattice dual to the one Yang-Mills lattice links are defined on^, the vortex world- 
surfaces are composed of elementary squares, and can thus link with the Wilson loops of the 
original lattice. As described and motivated in detail in |lllll2j . since the model is not intended 
to describe the structure of the theory at arbitarily short distances, it has a fixed lattice spacing 
a representing the vortex thickness. Continuity of the magnetic flux forces the vortex world- 
surfaces to be closed, and the statistical weight of a (closed) vortex world-surface is determined 
by a model action inspired by a gradient expansion of the Yang-Mills action of a center vortex 



Here, the (dual) lattice elementary square extending from the dual lattice site x into the positive 
IJ, and v directions is associated with a triality q^v{x) € {—1,0, 1} (in the case of an underlying 
SU{3) gauge group), labeling the center flux carried by that elementary square. The value 
Q^iui^) = indicates that the square is not part of a vortex world-surface, whereas the other 
two values are associated with the two types of center flux a vortex surface can carry^. The 
two terms in eq. (jlj correspond to a Nambu-Goto and a curvature term, respectively. They 
implement the notion that vortices, on the one hand, may be associated with a surface tension 
such that it costs a certain action increment e to add an elementary square to the surface. On 
the other hand, vortices are stiff, such that an action increment c is incurred for each pair of 
elementary squares in the vortex surface which share a link but do not lie in the same plane 

''Two lattices of spacing a are dual to one another if one can be generated by shifting the other one by the 
vector (1, 1, 1, l)a/2. 

■'To be precise, in terms of elementary Wilson loops (plaquettes) on the original lattice, the plaquette Uapiy) 
extending from y into the positive a and /? directions takes the value Uaf){y) = exp(i7r/3 • taPij.vqfi.vix)), where 
X = y + {ec + G/s ~ Bfj. — ev)a/2, with ex denoting the unit vector in the A-direction. 
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(i.e., vortex curvature is penalised). For details on the physical foundations of the model, and 
the determination of the couplings e and c, the reader is referred to refs. 112j . In practice, 
physical ensembles which reproduce the gross features of the corresponding Yang-Mills theory 
in the infrared can be achieved with e = 0. In the case of an underlying SU{2) gauge group, 
the value c = 0.24 generates a physical ensemble, whereas in the case of the SU{3) gauge group, 
the appropriate value is c = 0.21. These values were determined by requiring the models to 
correctly reproduce the ratio of the deconfinement temperature Tc to the square root of the 
zero-temperature string tension y/a found in the corresponding Yang-Mills theory. 
Since the random vortex world-surface model was proposed as an effective description of the 
long-range properties of Yang-Mills theory, including its deconfinement phase transition, the 
temperature dependence of &{T) near the transition is a non-trivial and important test of the 
model. In the random vortex world-surface model, the elementary squares q^ui-c) of the dual 
lattice represent the building blocks for the dynamical center vortex degrees of freedom; in this 
framework, it is thus particularly simple to introduce an additional vortex world-surface into 
the configurations. One simply needs to replace^ 

Qfiuix) — > {%v{x) + q) mod 3 (5) 

for all elementary squares q^u{x) making up the world-surface in question, with a fixed value of 
q G {—1, 1}. The two possible values of q are related by a space-time inversion which reverses 
the direction of magnetic flux (and leaves the action invariant); one can therefore restrict oneself 
to ^ = 1 without loss of generality. 

The dual string tension is obtained in this work specifically by calculating the excess free energy 
in the presence of a vortex world-surface P occupying an entire lattice plane extending into two 
spatial directions. The plane P is composed of elementary squares pi with !<«</, where / is 
the number of elementary squares in P. Let Q denote collectively a vortex world-surface con- 
figuration {q^y{x)}, and S[Q] the corresponding action @. It will be useful for the following to 
introduce a notation for introducing additional elementary vortex squares pi into configurations; 
this is best done recursively: Starting with a configuration Qq from the conventional random 
vortex world-surface ensemble, without any additional vortex surfaces introduced, define the 
configuration Qi as the configuration obtained by effecting the transformation © specifically on 
the elementary square pi of the configuration Qi^i. One can also define corresponding partition 
functions Zi as 

Zi = j [dQi\ exp(-5[Q,]) = j [dQo] exp(-5[Q,]) . (6) 

In the notation of eq. Zp = Z/. To obtain the dual string tension, it is thus necessary to 
evaluate 

Zi /[dQo]exp(-5'[Qo])exp(-(5'[Qj] - S[Qq])) _ 

-Tf- = rurTi / cm n = \'^^Pi~('^ V/ - ^[Qo )))o ■ ('j 

Zo J[dQo\e^p{-S[Qo\) 

However, this expression is of little practical use as it suffers from a serious overlap problem: The 
quantity being averaged varies over many orders of magnitude, so that most configurations only 
give an exponentially small contribution to the average. The resulting numerical noise precludes 
extracting a useful signal. This problem is addressed by using an algorithm introduced by 
de For cr and et al J7]. By decomposing 

Zi Zj Zi 



(8) 



^The modulo operation is to be applied such that the result again takes a value in { — 1, 0, 1}. 
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the problem is separated into the calculation of a (sizeable) number of independent expectation 
values with good overlap; namely, one evaluates the effect of introducing just one additional 
elementary square pi into the configurations at a time, 

_ J[dQo]eM-S[Q.-i])eM-{Sm - g[Q.-i])) ^ ^exp(-(g[Q.] - . (9) 



/[dQo] exp(-5'[Qi_i]) 

Note that the classes of configurations Qi, excepting Qq and Qj, contain open vortex world- 
surfaces, bounded by center monopoles, which violate the Bianchi identity. As already discussed 
further above, these are unphysical objects which are introduced here merely as intermediate 
mathematical constructs. If one implemented the Bianchi constraint in terms of an infinite 
action term penalizing center monopoles, all intermediate partition functions Zi (apart from Zq 
and Zj) strictly speaking would vanish. These Zi thus really are calculated with a modified 
weight in which the infinite self-energies of the center monopoles introduced explicitly into the 
configurations are left out (all other center monopoles are of course still completely suppressed; 
the integration /[(iQo] is over the conventional ensemble which respects the Bianchi identity). 
Of course, the final result for the dual string tension ultimately depends only on the physical 
partition functions Zj and Zq; all intermediate partition functions cancel out in the product Q. 
In practice, it is possible to improve the measurement of the individual expectation values @ 
further by using a multi-hit procedure. The quantity S[Qi] — S[Qi-i] depends only on the 
configuration in the neighborhood of the elementary square pi; thus, by performing multiple 
configuration updates and measurements in that neighborhood before carrying out the next 
global update of the configuration, statistics can be improved considerably without spoiling 
detailed balance or ergodicity. 

As the plane P is filled up one elementary square pi at a time^ , one can record the free energies 

Fi = -ln{Zi/Zo) (10) 

of the partial vortex surfaces, up to the free energy Fj of the full plane P. The dual string 
tension is extracted as 

*(T) = A (n, 

where a denotes the lattice spacing. As already mentioned above, the lattice spacing in the 
random vortex world-surface model is a fixed physical quantity related to the transverse thickness 
of the vortices; for a detailed discussion, cf. At the physical point e = 0, c = 0.21 of the 

SU{3) model, its value is a = 0.39 fm, where the scale was fixed by equating the zero-temperature 
string tension with (440 MeV)^. Besides the dual string tension, which is the principal quantity 
of interest in this work, below also the partial free energies Fi will be discussed. 



3 Numerical results 
3.1 General remarks 

Since the lattice spacing a of the random vortex world-surface model is taken to be a fixed 
quantity that is not scaled towards a continuum limit, the temperature in the model can only 
be changed in rather large steps if the coupling constants are kept fixed. However, by employing 
an interpolation method, a continuous range of temperatures can be explored |111 112j : For all 

^In practice, this was done row by row, as one would read a text. 
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measurements in this work, e = is fixed and only variations in the curvature coupling c are 
considered. For several temporal extensions A^o = 1)2,3 of the lattice, one can vary the coupling 
c until the deconfinement phase transition is observed at critical couplings c* = c*{Nq = i). The 
critical couplings for the gauge group SU (3) on lattices of extension Nq x 30'^ are^ 



No = i 


1 


2 


3 


c* 


0.0872 


0.2359 


0.335 



From the critical temperature aTc = 1/Nq for the three values c = c*, the function aTc{c) can be 
determined for all relevant couplings c by interpolation. This, in turn, specifies the temperature 
in units of the deconfinement transition temperature for any c and Nq via 

T _ 1 

¥c ~ No ■ aTc{c) ■ 

In particular, this permits carrying out measurements at a given T/Tc for different A^o and the 
corresponding c. Combining these measurements, the quantity in question, at the given T/Tc, 
can then also be obtained for the physical value c = 0.21 by interpolation in c. This procedure 
is used to arrive at the last two columns of Table ^ below. 

To be consistent with this determination of the temperature, and also in order to minimize 
finite-size effects, the evaluation of the dual string tension in this work was likewise carried out 
on A'^o ^ 30"^ lattices. To check for finite-size effects, control measurements were performed on 
lattices with spatial extensions of 10, 12 and 16 lattice spacings; in all cases, the discrepancies 
between measurements of the dual string tension at spatial extensions 16 and 30, respectively, 
were smaller than the statistical errors of the measurements. 

3.2 Boundary effects on open vortex world-surfaces 

The left-hand panel of Fig.^shows a typical result for the free energies Fi of partial vortex world- 
surfaces in the deconfined phase, at T/Tc = 1.093. The area of the surfaces in units of is i, 
and the approximate linear rise with i thus indicates an area law, as expected in the deconfined 
phase. However, there are also systematic deviations from a behavior purely proportional to the 
area, in particular for small i, and for large i, when the surface almost completely covers the 
lattice plane it is located in. 

These deviations can be formally attributed to self-energies and residual interactions of the cen- 
ter monopole world-lines which necessarily bound the intermediate open vortex world-surfaces. 
They are the most transparent when considering surfaces made up of full rows of elementary 
squares, i.e., when i is an integer multiple of the spatial extension of the lattice, i = kNg. In 
this case, there are two parallel center monopole world-lines winding along a spatial dimension 
of the lattice at a distance k ■ a and the deviations from the area law formally correspond to 
self-energies and residual interactions of the two monopoles. The right-hand panel in Fig. ^ 
displays the deviation AF{k) = -FfeAr^ — Fj ■ kNg/I from an area law as a function of the center 
monopole separation in lattice units k. It should again be noted that this quantity does not have 
an immediate physical interpretation. Apart from the monopoles propagating into a spatial in- 
stead of the temporal direction, they strictly speaking are associated with an infinite self-energy 

^The value of C3 was only determined to within an error of 1% due to its weak influence on tfie interpolations 
needed in practice. The other two values are accurate to the digits shown. 
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Figure 1: Left panel: Free energies Fi, cf. eq. (|lfl|) . of partial vortex world-surfaces, as a function 
of the number i of elementary squares included in the surface. Measurements were taken on 
a 30'^ X 2 lattice at e = and c = 0.2550, which corresponds to T/Tc = 1.093 (deconfined 
phase). Right panel: The residual magnetic monopole correlation AF{k) = i^fcAr^ — Fj ■ kNg/I 
as a function of the monopole pair separation in lattice units k, from the same measurement. 
The curve represents a fit by a (self-energy) constant plus a Yukawa potential, AF{k) = Vq + 
exp{—fj.k)/k. 



preventing their appearance in the physical vortex ensemble. An additional finite contribution 
to this infinite self-energy is physically inconsequential. 



3.3 The dual string tension 

Approaching the deconfinement phase transition temperature Tc from below, one can verify that 
the dual string tension is compatible with zero. Corresponding measurements on 30^ x 2 lattices 
at couplings c below the critical coupling c = 0.2359, which realize the confining phase, yield: 



c 




0.21 


0.000190 ± 0.000214 


0.23 


3.6- 10~^ ±0.001 


0.235 


-0.00037 ± 0.001 



Also, for comparison, a measurement on a 16'^ x 2 lattice at c = 0.21 resulted in aa^ = 
-0.000205 ± 0.000271. 

Above the phase transition temperature, measurements were carried out for a sequence of tem- 
peratures T/Tc. As displayed in Tabled for each temperature, measurements were performed 
both on a 30^ x 2 and a 30^ x 1 lattice with the appropriate couplings c realizing the given 
temperature. Thus, for a fixed temperature T/Tc, the dual string tension a for the two values 
ci {Nq = 1) and C2 (A'^o = 2) is obtained. To arrive at the physical value, those two results are 
then interpolated linearly to the physical coupling Cphys = 0.21. Since the lattice spacing at the 
physical point is known in absolute units, the physical value of the dual string tension can also 
be given in absolute units. This is shown in the last column of Table ^ These results are also 
presented graphically in Fig. El 
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No = l 


0.1073 


188 


134 


195 


No = 2 


0.2550 


110 


1.118 


No = l 


0.1123 


226 


164 


216 


No = 2 


0.2600 


133 


1.326 


No = 1 


0.1495 


437 


385 


330 


iVo = 2 


0.3000 


308 


1.500 


No = l 


0.1756 


513 


486 


371 


No = 2 


0.3350 


387 


1.629 


No = l 


0.1929 


548 


540 


391 


No = 2 


0.4500 


435 


1.770 


No = l 


0.2100 


579 


579 


405 


No = 2 







Table 1: Dual string tension a measured on 30"^ x A^o lattices with e = 0. The curvature couplings c 
were chosen such that the corresponding measurements for iVo = 1 and A^o = 2 implement the same 
temperature T/Tc- The results in the last two columns represent the interpolation to the physical value 
c = 0.21. Since the lattice spacing at this point is fixed to a = 0.39 fm, the value of the physical dual 
string tension (iphys can be given in absolute numbers. 
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Figure 2: Dual string tension a as a function of the temperature T in the deconfined phase, 
taken from Tabled 



As can be seen from the data, the dual string tension a rises with increasing temperature 
in the deconfined phase. As the phase transition is approached from above, the dual string 
tension a quickly vanishes as is expected for an order parameter. Since the phase transition 
in the SU{3) case is weakly first order |12j . a weak discontinuity is expected in the dual string 
tension <t(T). The discontinuity can be obtained by using the two data points nearest to Tc and 
extrapolating linearly to T = Tc. The statistical uncertainty for the dimensionless combination 
^■phys (30a)^ (the second to last column in Tab. ^ is 6.3 it 0.4 at T/Tc = 1.0015 and 11.5 it 0.5 
at T/Tc = 1.0052. Varying the slope of the extrapolation line within these error bounds, the 
resulting value at T = Tc+ is 



phys 



T=Tc + 



(34.5 ± 3.3) MeV , 



(12) 



As a cross-check, a spline extrapolation through the five data points nearest to the transition 
(up to T/Tc = 1-0195) yields ^^CTphys I t=Tc+ ^ 36.8 MeV, whereas using only the three nearest 
data points leads to ^/^phys I t=Tc+ ^ 36.5 MeV, both consistent with ((T^ . Also, linear regres- 
sion through the nearest three data points yields \/crphys|T=rc+ ~ 35.5 MeV. The data thus 
provide a clear signal for the discontinuity of the dual string tension a at T = Tc, in accor- 
dance with the first order character of the phase transition previously verified for the SU{3) 
random vortex world-surface model via the action density distribution at the critical tempera- 
ture . The small value of the discontinuity compared with typical strong interaction scales 
again demonstrates the weak first order character of the transition. 



4 Discussion 

The dual string tension in the framework of the SU{3) random vortex world-surface model was 
verified in this investigation to represent an order parameter for the deconfinement phase transi- 
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tion, and to furthermore reflect the weak first order character of the transition. Quantitatively, 
the result ()12[) for the discontinuity at the transition, 



= (34.5 ± 3.3) MeV , 

is smaller than one would expect from measurements of the order-disorder interface tension in 
SU (3) lattice Yang-Mills theory assuming perfect wetting. The available data for the order- 
disorder tension [2^ indicate a value of cxod ~ 0.015 T^, which translates to aod ^ (34MeV)^, 
if the scale is set by equating the zero-temperature string tension with (440 MeV) ^. Assuming 
perfect wetting, the order-order interface tension, equivalent to the dual string tension, would 
be expected to take the value 2ao(i w (48MeV)^. 

There are, however, caveats to this comparison. On the one hand, though wetting has been 
observed in SU{3) lattice Yang-Mills theory [261 I27j . there remains an uncertainty as to its 
quantitative extent. On the other hand, there are also still significant uncertainties involved in 
extrapolating the aforementioned measurements of the order-disorder tension from the lattice 
spacings they were obtained at to the continuum limit. Current efforts to accurately determine 
the dual string tension at the deconfinement transition in SU (3) lattice Yang-Mills theory j21j 
actually indicate significantly higher values of the tension at coarse lattice spacings; however, 
also in this case, no conclusive statement concerning the extrapolation to the continuum limit 
can be made on the basis of the preliminary data presently available. These data display a 
substantial downward trend for the dual string tension as the lattice spacing is reduced. 
In view of the current status, it seems premature to draw a final conclusion as to whether the 
result obtained here within the random vortex world-surface model quantitatively reflects the 
behavior of full SU{3) Yang-Mills theory or not. On the other hand, it should be noted that it is 
possible to significantly enhance the first-order character of the deconfinement phase transition 
in the vortex model by introducing an additional term into the action which favors vortex 
branchings; the latter are presumably instrumental in establishing first-order behavior, since 
they represent the feature which qualitatively distinguishes SU{3) vortex configurations from 
SU{2) vortex configurations. Introducing such a mechanism has indeed proven to be efficacious 
in preliminary investigations of the S'C/(4) random vortex world-surface model, which will be 
reported in a separate publication. Whether it will ultimately be necessary to invoke this type of 
mechanism in the SU (3) case is as yet unclear in view of the current status of the data discussed 
above; however, in the S'C/(4) case, such an additional term in the action will certainly play a 
role. 
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